CrownDiameter Function

public function CrownDiameter(dbh, den, denmin, denmax, dbhdcmin, dbhdcmax) result(cd)

compute crown doameter (m)

References:

Collalti, Alessio & Perugini, Lucia & Santini, Monia & Chiti, Tommaso & Nolè, Angelo & Matteucci, Giorgio & Valentini, Riccardo, 2014. A process-based model to simulate growth in forests with complex structure: Evaluation and use of 3D-CMCC Forest Ecosystem Model in a deciduous forest in Central Italy, Ecological Modelling, 272(C), 362-378.

Arguments

Type IntentOptional Attributes Name
real(kind=float), intent(in) :: dbh

diameter at brest height (cm)
real(kind=float), intent(in) :: den

current tree density (trees/ha)
real(kind=float), intent(in) :: denmin

minimum tree density (trees/ha)
real(kind=float), intent(in) :: denmax

minimum tree density (trees/ha)
real(kind=float), intent(in) :: dbhdcmin

minimum ratio between stem and crown diameters (m/cm)
real(kind=float), intent(in) :: dbhdcmax

maximum ratio between stem and crown diameters (m/cm)

Return Value real(kind=float)

returned result


Variables

Type Visibility Attributes Name Initial
real(kind=float), public :: dbhdc

actual ratio between dbh and crown diameters (m/cm)

Source Code

FUNCTION  CrownDiameter &
!
(dbh, den, denmin, denmax, dbhdcmin, dbhdcmax) &
!
RESULT (cd)


IMPLICIT NONE

!arguments with intent(in):
REAL (KIND = float), INTENT(IN) :: dbh !! diameter at brest height (cm)
REAL (KIND = float), INTENT(IN) :: den !! current tree density (trees/ha)
REAL (KIND = float), INTENT(IN) :: denmin !! minimum tree density (trees/ha)
REAL (KIND = float), INTENT(IN) :: denmax !! minimum tree density (trees/ha)
REAL (KIND = float), INTENT(IN) :: dbhdcmin !! minimum ratio between stem and crown diameters (m/cm)
REAL (KIND = float), INTENT(IN) :: dbhdcmax !! maximum ratio between stem and crown diameters (m/cm)

!local declarations:
REAL (KIND = float) :: cd !! returned result
REAL (KIND = float) :: dbhdc !! actual ratio between dbh and crown diameters (m/cm)

!-----------------------------------end of declarations------------------------

!compute actual dbhdc
dbhdc = DBHDCeff (dbhdcmin, dbhdcmax, den, denmin, denmax)

cd = dbh * dbhdc

RETURN
END FUNCTION CrownDiameter